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1 Introduction 


There are many engineering problems that involve linear systems of equali¬ 
ties and inequalities. Engineers are familiar with systems of equalities with a 
unique solution, but many of them are not used to deal with systems of equa¬ 
tions with multiple solutions and their interpretations. The problem is even 
worse when systems of inequalities are dealt with, because only few people 
know how to obtain the associated compatibility conditions and solve them. 
These two types of systems can be interpreted from the mathematical or the 
engineering points of view, which are complementary and provide a deep un¬ 
derstanding of the problem under study. However, people working in these 
problems use to have knowledge about only one of these two perspectives and 
are unaware of the relations between the mathematical and the engineering 
concepts, which leads to important limitations in the capacity of extracting 
conclusions from the results that can be expected after a careful analysis of 
these problems from both points of view. 


This paper points out these relations and makes them explicit for the readers 
to discover the new world that arises when contemplating the compatibility 
conditions or the set of general solutions from this dual perspective. 


In this paper, we have selected a particular example to illustrate these two 
points of view, the truss analysis problem, and we exploit this dual (math¬ 
ematical and engineering) perspective to deal with a problem that involves 
linear systems of equalities or inequalities, depending on the constraints used 
to model the reality. As we shall see, many questions of practical interest arise 
and can be answered thanks to this dual analysis of the problem. 


The paper is structured as follows. In Section 2 the problem of determining the 
compatibility conditions of systems of linear equalities and solving them to¬ 
gether with an analysis of the general mathematical structure of their solutions 
is discussed. In Section 3 the truss analysis problem is described and used to 
illustrate all the theoretical methods. In Section 4 we classify truss structures 
according to some mathematical criteria in isostatic, hyperstatic, mechanism 
and critical trusses. In Section 5 we discuss the compatibility of systems of lin¬ 
ear inequalities and the mathematical structure of their solutions. In Section 
6 we revisit the truss analysis problem adding some constraints to illustrate 
the methodology and several engineering problems. Finally, in Section 7 some 
conclusions and recommendations are given. 
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2 Dealing with Systems of Equations 


In many engineering applications we find systems of linear equations of the 
form: 


Ax = b, (1) 

where A G IR' nxn , and b G IR m , being m and n the number of equations 
and unknowns, respectively, which can be equal or not. Before trying to solve 
a system of equations of this type it is interesting to check whether or not 
the system is compatible, i. e., if it has some solutions. These conditions can 
be given in terms of bi, b 2 ,..., b m , and always have a physical or engineering 
meaning that raises some light about the problem under consideration. 

Once the system is proven to have solution we can obtain the set of all its 
possible solutions. To solve these two problems an algorithm that gives the 
linear space orthogonal to the linear space generated by a set of vectors can 
be used (for a detailed description of this algorithm see Castillo, Cobo, Ju- 
bete, Pruneda and Castillo [3], Castillo, Cobo, Fernandez-Canteli, Jubete and 
Pruneda [1], Cobo, Jubete and Pruneda [2]). 


2.1 Deciding whether or not a linear system of equations is compatible 


The system (1) can be written as 


aqai + x 2 a- 2 H-b x n a n = b, (2) 

which shows that the vector b = (b\b m ) T belongs the linear space gener¬ 
ated by the column vectors {a!, ai, ..., a ra } of the system matrix A, i.e., the 
compatibility requires: 


b G £(ai, a 2 ,..., a n ) ^bG ^£(ai, a 2 ,..., an) -1 ^ , 

where J_ refers to the orthogonal set. Thus, analyzing the compatibility of the 
system of equations (1) reduces to finding the linear subspace C (wi,..., w p } 
orthogonal to C {ai,..., a re } and checking whether or not h 1 W = 0. 

Example 1 (Compatibility of a linear system of equations) Suppose that 
we are interested in determining the conditions under which the system of 
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equations 


x x +x 2 -x 3 +x 4 +x 5 = a 
—x 2 +2x 3 +x 4 -2x 5 = b 

X\ —x 2 +2x 3 = c (3) 

— 2xi +2x 2 —3x 3 +x 4 —x§ = d 
x\ +x 2 +x 3 +x 4 +x 5 = e 

is compatible. Then, first we obtain the linear subspace orthogonal to the 
linear subspace generated by the column vectors in (3) that is: 

W = £ { Wl } = £ {(0, 2, -7, -3,1) T } , (4) 

which implies the following compatibility condition: 

wf (a, b, c, d, e) T = (0, 2, —7, —3, l)(a, b, c, d, e) 1 — 0 =>• 2b — 7c — 3d + e = 0. 

( 5 ) 


2.2 Solving a homogeneous system of linear equations 


Consider the homogeneous system of equations 

Ax = 0, 


which can be written as 


a*x 7 =0; i = 1,..., m. 


( 6 ) 

(V 


Expression (7) shows that (x 4 ,, x n ) is orthogonal to the set of row vectors 

{a 1 , a 2 ,..., a m } of A. 

Then, obtaining the solution to system (7) reduces to determining the linear 
subspace orthogonal to the linear subspace generated by the rows of matrix 
A. Thus, the general solution of an homogeneous system of linear equations 
is a linear space, i.e., of the form 

p 

x = Pi e 1R. 

i =1 

Example 2 (An homogeneous system of linear equations) The system 
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of linear equations 


X\ +x 3 — 2x 4 +x 5 = 0 

x 1 -x 2 -X 3 +2x 4 +X 5 = 0 

—X\ +2x 2 +2x 3 —3x4 +x 5 = 0 


has as general solution the linear space 




M 


f°) 

^2 


2 


2 

x 3 

= Pi 

-3 

+ P2 

-7 

X 4 


-1 


-3 

K x 5/ 


v 0 J 


v x y 


Pi, P-2 £ IR ¬ 


IS) 


2.3 Solving a complete system of linear equations 


Now consider again the complete system of linear equations (1), that adding 
the artificial variable x n+ 1 , can be written as 


+ 0 \ 2 x 2 + ■ ■ ■ +ai n x n —bix n +i — 0, 

Q' 2 i^i +a 22 x 2 + ■ ■ ■ +a 2n x n —b 2 x n +i = 0, 

+Q J m2X 2 + +a mn X n b m X n -|_i 0, 

Xn +1 1 - 


(9) 


The first m equations of the system (9) can be written as 


($11, • • • 5 ni i *^n+l) b? 

(&21> • • • j ^2 n't ^2)(«^1 ? • • • j ^n+l) 0, 


(flml) • • • 5 • • • ? *£715 ^n+l) 0 ? 


( 10 ) 


which shows that (aq,..., x n , x n+1 ) is orthogonal to the set of vectors 

\(®11 j • • • j ®1 n, ^1 ), (®21 , • • • , ®2 n, ^ 2 ), • • • , (®ral , • • • , ®mn? ^m)/j 
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i.e., the solution of (10) is the linear subspace orthogonal to the linear subspace 
generated by the rows of matrix A b , i.e., 

£{(an, . . . , ai n , —bi), ((Z2Ij • • • , «2 ni ^2), • • • , («ml )•••■) ^mni ^m) } 


Thus, the solution of (1) is the projection on Xi x • • • x X n of the intersection 
of such a orthogonal linear subspace and the set {x|x n+1 = 1}. In other words, 
the general solution of a complete system of linear equations is an affine space, 
that is the sum of a constant vector plus a linear space: 

p 

x = x° + ^2 Pi v it Pi e H, 

i= 1 

where the first vector in the right hand side is an arbitrary particular solution 
of the system (1), and the linear space is the set of solutions of the associated 
homogeneous system. 


Example 3 (A complete system of linear equations) The system of lin¬ 
ear equations 

x 1 +x 3 -2 x 4 +x 5 = 1 

x\ -x 2 —x 3 +2 x 4 +£5 = 2 (11) 

—x\ +2^2 +2x 3 —3x4 +X5 = 1 
has as general solution the affine space 



Pi , P 2 £ 1R • 


3 The Truss Analysis Problem 


Consider the engineering problem of structures consisting of trusses, i.e., struc¬ 
tures made of bars joined by frictionless hinges which imply that only tension 
or compression stresses exist in the bars. A typical truss structure is shown in 
Figure 1. 

With reference to the bar element (see Figure 2) that has length L and cross- 
section A, let us first establish the stiffness k. By definition, the force N is 
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Fig. 1. Two-dimensional truss structure consisting of bars joined by frictionless 
hinges. 


N 


© 

0 = 


© 

=0- 


x' 


Fig. 2. Loads acting on an individual bar. 
related to the normal stress a as 

N = Aa, (12) 

and assuming that the bar behavior is linearly elastic, according to Hooke’s 
law, we have 

a = Ee , (13) 

where E is the Young’s modulus and e is the longitudinal strain, that is as¬ 
sumed constant along its length, and has value 


where u! i and «'• are the longitudinal displacements of the end points i and j, 
respectively. 

Combining Equations (12)-(14) we obtain 

N = k(uj — u[), (15) 

where k = AE/L is the bar stiffness. 



To calculate a given truss structure first we need to know the behavior of 
a single bar in an arbitrary direction. Thus, consider the bar shown in Fig¬ 
ure 3 (a). Our purpose is to determine the relation between the node forces 
F Xi , F yi , F Xj , F yj and the node displacements Ui,Vi, Uj, vj. 

To this aim, we introduce a local x'-y'-coordinate system (see Figure 3 (b)) 
with its origin at the end point i and its axis x' directed along the bar axis 
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Fig. 3. Displacements and force components of a bar element with respect to: (a) 
the global x-y system, and (b) the local x'-y' system. 

from the end point i towards the end point j. Note that this coordinate system 
can be obtained from the ^-system by a rotation of angle a. 


From Figures 2 and 3 (b) and Equation (15) we obtain 


K t = -N = k(u[ - «') 
K, = N = k(u' t - «'), 

which can be written as 


( P' \ 

Xi 


' 1 0-10^ 



F' 


0 0 0 0 


V 'i 

Vi 

= K e u e = k 




F' 

Xj 


-10 10 





o 

o 

o 

o 




(16) 


(17) 


where YF' is the so-called local element stiffness matrix, that is symmetric. 
It can be observed that the displacements along the r/'-axis, i.e., v- and v'j, 
produce no forces on the bar. This assumption is valid only if the displacements 
are small; otherwise, the components v[ and Vj might result in an elongation 
of the bar and thus in the development of forces. 

Let us now establish a relation between the displacements Ui,Vi,Uj,Vj and 
?j(, v' i: u'p v'y By geometrical arguments it follows directly that 


Ui = u[ cos a — v[ sin a 
Vi = u[ sin a + v\ cos a 


( 18 ) 


and the corresponding equations for u 3 and Vj. 






















Therefore 




' COS CK 

— sin a 

0 

0 ' 



Vi 

T / 

sin a 

cos a 

0 

0 


v\ 


= 1 / u e = 







Uj 


0 

0 

COS OL 

— sin a 


u 'j 



l 0 

0 

sin a 

cos a j 


{ ) 


( 19 ) 


where L e is the so-called transformation matrix, that is an orthogonal matrix 
(L e_1 = I/ T ). 


Premultiplying (19) by L e and using the orthogonal property, we then conclude 


that 


u e = L e u e 


( 20 ) 


As both forces and displacements are vector quantities, we can choose the 
same basis for the forces and then, the element forces F e and F e are related 
in exactly the same manner as the element displacements, i.e., 

F e = L ,;T F /e . (21) 


With these preliminary remarks, we can determine the relation between the 
element forces F e and the element displacements u e . Insertion of (20) into (17) 
yields 

F e ' = K e L e u e , (22) 

and premultiplication by L e ' and using (21) leads to 

F e = L eT F ,e = L e K c L e u e = K e u e , (23) 


where K e is the so-called global element stiffness matrix. Its explicit expression 

is 


K e 


( 2 • 2 • \ 
cos a cos asm a — cos a — cos a sin. a 

cos a sin a sin 2 a — cos a sin a — sin 2 a 

— cos 2 a — cos a sin a cos 2 a cos a sin a 

^ — cos a sin a — sin 2 a cos a sin a sin 2 a j 


(24) 


which shows the symmetric character of the global element stiffness matrix 

K e . 


The truss analysis problem has the following elements to be considered: 

Bars: The longitudinal elements supporting only tension or compression 
stresses (the number of bars is b). 
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Nodes: The frictionless elements that join the bars allowing relative rotation 
between them (the number of hinges or nodes is m). 

Supports: Elements that prevent the structure to experiment solid rigid 
displacements. The reactions exerted by the supports are the necessary ones 
for the equilibrium to hold. The number of constraints associated with the 
supports is denoted c. 

Forces: The forces acting on the nodes (hinges) of the structure that are the 
data of our problem. 

Displacements: Movements suffered by the nodes owing to the strains pro¬ 
duced by the stresses in the bars. 

Unknowns: The hinge (nodal) displacements. The number of unknowns co¬ 
incides with the number of possible hinge displacements, which is the num¬ 
ber of degrees of freedom (two times the number of nodes m in case of 
bi-dimensional analysis). 

Equations: The mathematical relations that give the nodal forces in terms 
of the nodal displacements. To derive the system of equations that model a 
given problem the global structure stiffness matrix must be obtained. 



Example 4 (A simple truss structure) As an example of how the ele¬ 
ment stiffness matrix can be used to derive the stiffness matrix for the whole 
structure, consider the simple truss structure shown in Figure 1. Using the 
relationships between forces and movements, establishing the equilibrium of 
vertical and horizontal forces for each node, and replacing the static magni¬ 
tudes (stresses) by their equivalents as a function of the nodal displacements, 
then the so-called global assembled stiffness matrix is obtained. Note that 
each element stiffness matrix (K e ) is assembled into the global element stiff- 
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ness matrix (K), to get the following system of equations: 


Ku = 


150 


' 179 

72 | - 

-125 

0 | 

-54 

-72 | 

0 

0 

72 

96 | 

0 

0 

-72 

-96 | 

0 

0 

— 

- + 

— 

- + 

— 

- + 

— 

— 

-125 

0 | 

233 

o 1 

-54 

72 | 

-54 ■ 

-72 

0 

0 

0 

192 | 

72 

-96 

-72 - 

-96 

— 

- + 

— 

- + 

— 

- + 

— 

— 

-54 - 

-72 | 

-54 

72 | 

233 

0 | - 

-125 

0 

-72 - 

-96 

72 

-96 

0 

192 

0 

0 

— 

- + 

— 

- + 

— 

- + 

— 

— 

0 

o 1 

-54 

-72 | - 

-125 

o 1 

179 

72 

0 

o 

-72 

-96 

0 

o 

72 

96 


Ml 


F 

± XI 

V\ 


F 

1 yi 

u 2 


F 

A X2 

V2 


r S 1 

^ 1 

u 3 


F 

± X3 

V-3 


F 

1 V3 



f X4 

, VA ) 


k F y 4 j 


= F 


(25) 


where Ui, % = 1,..., 4 are the node displacements shown in Figure 4 (a) and 
F Xi , Fy t ; i = 1,..., 4 are the external forces acting on the hinges of the structure 
including support reactions (see Figure 4 (b)). The different partitions refer 
to different nodes. 


First, using the method described in Section 2.1 the conditions to be satisfied 
by the linear system (25) to have solution (one or more) are obtained: 


0 — F X1 + F X2 + F X3 + F Xi (26) 

0 = Fyi + Fy* + Fyz + Fy* ( 27 ) 

0 = 4/5 F X1 - 3/5 F yi + 4/5 F X2 + 3/5 F y2 + 6/5 F y4 , (28) 

where (26) and (27) express the equilibrium of horizontal and vertical forces 
respectively, and (28) establishes the equilibrium of moments (the moment are 
taken with respect to node 3). Thus, these conditions state that the structure 
subject to the external forces and the support reactions must be in equilibrium 
(E F x . =0,E4 = 0 and E = 0). 

The rank of K in (25) is 5 that implies that three (8 — 5) rows are linear 
combinations of the other rows, and then the system (25) has infinite solutions. 
So, some extra equations of the form Bu = b, where b is the boundary 
conditions vector, can be added maintaining the compatibility of system (25). 

To obtain the set of all possible solutions of system (25), we apply the orthog- 
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onalization algorithm and obtain: 


H 


1 467/120 N 


M 


fo) 


1 4/5 N 

Vi 


501/160 


0 


1 


-3/5 

U 2 


359/120 


1 


0 


4/5 

V2 

1 

-359/160 


0 


i 


3/5 


7 . 

+pi 


+ P2 


+ P3 

U'i 

k 

39/10 


1 


0 


0 

V 3 


0 


0 


1 


0 



0 


1 


0 


0 

UJ 


v 0 y 


v°y 


v 1 2 3 * ; 


v 6 / 5 / 


Pi) Pli P3 £ 1R, 


( 29 ) 

which from a mathematical point of view is an affine space, i.e., the sum of 
a given vector (the first one) plus a linear space of dimension 3 (an arbitrary 
linear combination of 3 linearly independent vectors). 



Fig. 5. Illustration of the basis of the linear space of dimension three that appears in 
the general solution of the two-dimensional truss structure in Example 5. The basis 
has three generators that correspond to: (a) a horizontal translation, (b) a vertical 
translation, and (c) a rotation. 

Note that the linear space is generated by three vectors that, due to the 
symmetry of the stiffness matrix, coincide with those in (26)-(28). From an 
engineering point of view, this solution must be interpreted as follows: 


(1) The dimension of the linear space is the number of degrees of freedom 
we have, i.e., the maximum number of u and v displacements that can 
be fixed for the system to have a unique solution. 

(2) The linear space component of the solution (unlimited values of the p 
coefficients, and consequently unlimited displacements) implies that the 
structure can be located anywhere whereas the forces acting on the struc¬ 
ture guarantee equilibrium. 

(3) The first vector is a particular solution, i.e., a solution to the stated 

problem. Note that it satisfies Equation (25) for the F values in Figure 6. 

It is worth noting that this vector can be replaced by any other particular 
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solution, which can be obtained by adding to it any linear combination 
of the three basic vectors. 

(4) The second vector corresponds to a solution of the associated homoge¬ 
neous problem, i.e., with no external forces. In this particular case it 
corresponds to the horizontal displacement of the structure as a rigid 
solid (see Figure 5(a)). 

(5) The third vector corresponds to another solution of the associated ho¬ 
mogeneous problem. In this case, it is the vertical displacement of the 
structure as a rigid solid (see Figure 5(b)). 

(6) Finally, the fourth vector corresponds to another solution of the asso¬ 
ciated homogeneous problem, which corresponds to the rotation of the 
structure as a rigid solid with respect node 3 (see Figure 5(c)). 

It is obvious that the linear space generated by the last three vectors in (29) 
can be represented using another basis of the same space, i.e., considering 
equilibrium with respect other point. 


4 Classification of truss structures 


Structures can be classified in different ways depending on the rank of the 
matrix ^ K | B ^ and the value of b + c and m, as shown in Table 1. They 
include isostatic, hyperstatic, critical and mechanism truss structures: 


Isostatic: In this kind of structures 2m = b + c and the rank ^ K | B 

2m, i.e., a unique solution exists. They are characterized because (a) the 
boundary conditions ensure equilibrium under all possible external forces, 
(b) thermic strains do not induce stresses, and (c) if any of their elements 
(nodes, bars, or boundary conditions) is removed the structure becomes a 
mechanism or a critical structure. 

Hyperstatic: In this kind of structures b + c > 2m where (b + c) — 2m is 
the degree of hyperstaticity and the rank | = 2m, i.e., a unique 

solution exists. They are characterized because (a) the boundary conditions 
ensure equilibrium under all possible external forces, (b) thermic strains 
induce stresses, and (c) if any of their elements (nodes, bars, or boundary 
conditions) is removed the structure becomes a isostatic structure or remains 
hyperstatic. Note that there are two hyperstatic cases, one is due to an 
excess of bars in the truss structure, and the other is due to an excess in 
the supports (boundary conditions). 

Mechanism: In this kind of structures b+c < 2m and the rank ( K I B ^ < 



2m, i.e., infinite solution exists. They are characterized because (a) the 
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Table 1 

A classification of truss structures. 



/ K \ 

Rank — 

V B / 

b + c 

Isostatic 

2m 

2m 

Hyperstatic 

2m 

> 2m 

Critical 

< 2m 

> 2m 

Mechanism 

< 2m 

< 2m 


boundary conditions do not ensure equilibrium under all possible external 
forces, (b) thermic strains do not induce stresses. 

Critical: In this kind of structures 2m > b+c and the rank ^ K | B ^ < 2m, 

i.e., infinite solution exists. They are characterized because (a) the boundary 
conditions do not ensure equilibrium under all possible external forces within 
the hypothesis of small deformations, (b) thermic strains induce stresses, 
and (c) if any of their elements (nodes, bars, or boundary conditions) is 
removed the structure becomes a mechanism. 

To illustrate their main differences we include some examples. 

Example 5 (An isostatic truss structure) Consider the particular case 
shown in Figure 6, where there are 3 external forces. Since the vector of exter¬ 
nal forces F has to fulfill the compatibility conditions (26)-(28), we use them 
to obtain the truss reactions: 



Fig. 6. The two-dimensional truss structure showing the external forces and reac¬ 
tions. 


F x i=9/4; F yi — 2; F Xi = -13/4; 
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T 


and then, the force vector becomes: F = 1/4 (^9 80 —4 4 —4 —13 OJ • 

Next, we consider the boundary conditions imposed by the supports that 
establish the final location of the structure. Assuming that those boundary 
conditions are ip = 0,i>i = 0 and u 4 = 0, then from (29) we have: 


u 1 = 467/(120/c) + pi + 4/5p 3 = 0 
V\ = 501/(160fc) + p 2 - 3/5p s = 0 
v 4 = p 2 + 6/5p 3 = 0. 


The solution of this system of equations is p\ = 0, p 2 = —121/(20 k) and 
p 3 = —467/(96A;). Thus the corresponding solution obtained after replacing 
these values of pi,p 2 and p 3 into (29), is unique and equal to: 


f«A 


0 1 

Vl 


0 

U 2 


- 9/10 

V2 

1 

-897/80 

u 3 

“ k 

39/10 

V-3 


-121/120 

U 4 


0 

K v *) 


v -951/80 j 


Note that the uniqueness of solution implies that the imposed boundary con¬ 
ditions are enough to avoid rigid solid movements and therefore the supports 
would be able to react to external forces acting on the structure in such a way 
that the compatibility conditions (equilibrium) hold. The solution is unique 
because for the submatrix corresponding to p coefficients and U\,v 4 and v 4 in 
(29), we have: 


Rank 


( l 0 4/5 ^ 
0 1 -3/5 
v 0 1 6/5 J 


= 3 




rank 




v B / 


2 m = 8. 


Example 6 (An isostatic truss structure) Consider the structure shown 
in Figure 7 subject to three forces Pi,P 2 and P 3 . Applying the superposition 
principle and assuming that the forces acting on the structure are: 
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Fi=Pi (o 0 0 0 1 0 -1 o) , 
f 2 = P 2 (3/4 1000-1-3/40) 

Fs = P 3 [ 3/2 10-100 -3/2 0 ) 

F 4 = P 4 (9/4 1 0 0 0 0 -9/4 -l) 

that satisfy the compatibility conditions (26)-(28), the associated system of 
equations leads to the general solution: 



( 3 ^ 


( 9 ^ 


/ 341 \ 


/287\ 



( 4 \ 

( u \ 


5 


20 


210 


60 


(A 


fo) 


5 

LL\ 


9 


19 


25 


27 





3 

Vi 


20 


To 


32 


^0 


0 


1 


5 


3 


9 


233 


179 






4 

u 2 


5 


20 


120 


60 


1 


0 


5 

V 2 

U 3 

Pi 

k 

9 

11 

k 

27 


233 

k 

27 


0 

1 


1 

0 


3 

20 

6 

80 

9 

Jr 3 

1 k 

160 

9 

40 

9 

+Pl 

+ P2 

+P3 

5 

0 



5 


10 


5 


5 






v 3 


0 


0 


0 


0 


0 


1 


0 

U4 










1 


0 


0 


0 


0 


0 


0 






W 


v° J 


l 0 J 


l 0 J 


v° J 


0 


1 


6 







' 5 / 


where Pi, P 2 , P 3 , Pa, Pi, Pi, p 3 £ K. From a mathematical point of view this 
is a linear space, i.e., the solutions are the linear combinations of seven given 
vectors. 



Fig. 7. Isostatic structure in Example 6 showing the applied loads. 

Example 7 (A critical truss structure) If we consider the same structure 
as in Example 5 with other boundary conditions, for example, U\ = 0, V\ = 0 
and u 2 = 0, then from (29) we have: 
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Ui = 467/(120A;) + p 1 + 4/5p 3 = 0 
Vl = 501/(160A:) + p 2 - 3/5p 3 = 0 
u 2 = 359/(1207c) + Pl + 4/5p 3 = 0, 


which is an incompatible system of equations, because for the submatrix as¬ 
sociated with the p’s and Ui,Vi and u 4 in (29) we have 


Rank 


'l 0 4/5 ^ 
0 1 -3/5 
V 1 0 4/5 y 


2 < 3. 


Adding the equations imposed by the boundary conditions to the system (25), 
considering the external forces (see Figure 8) F = 1/8 ^9 12 0 —8 8 —8 —17 4^ 

that satisfy the compatibility conditions (26)-(28) (equilibrium conditions) and 
applying the orthogonalization algorithm the following set of all possible so¬ 
lutions is obtained: 


f«A 


° 1 


0 \ 

Vl 


0 


0 

U 2 


0 


0 

V2 

1 

-49/16 


6/5 


= — 


+ P 


u 3 

k 

-1/24 


-4/5 

V-3 


-37/16 


3/5 

Ui 


-73/24 


-4/5 

K V,j 


V 0 J 


v 9/5 J 


(31) 


From an engineering point of view, this solution must be interpreted as follows: 

(1) The first vector is a particular solution, i.e., a solution to the stated 
problem. Note that it satisfies equation (25) for the F values considered 
and the boundary conditions u\ = 0,Ui = 0 and u 2 = 0. 

(2) The second vector corresponds to a solution of the associated homoge¬ 
neous problem, and represents the rotation of the structure as a rigid 
solid with respect node 1, as it is shown in Figure 8. This implies that 
there is a set of infinite solutions (rotations), and that the actual bound¬ 
ary conditions are not able to avoid the rigid solid displacement of the 
structure without the development of new forces in the bars under the hy¬ 
pothesis of small displacements. Moreover, the supports are not capable 
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Fig. 8. Illustration of the set of general solutions associated with the boundary 
conditions u\ = 0, v\ = 0 and U 2 = 0 for the critical two-dimensional truss example. 


of supplying the reactions required to satisfy the equilibrium conditions 
for all possible forces acting on the nodes of the structure. 



Fig. 9. The mechanism structure analyzed in Example 8. 


The same result is obtained considering 


B 


1 \ 0 0 0 0 0 0 o N 
01000000 
v 0 0 1 0 0 0 0 o y 


rank 




\ b j 


2m — 1 = 7 


that implies that there is one rigid solid rotation (p) allowed. 


Example 8 (A Mechanism structure) Consider the simple mechanism truss 
structure shown in Figure 9. The system of equations resulting from assern- 
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bling the global element stiffness matrices becomes: 


Ku 


' 54 72 0 0 -54 -72 0 0 

72 96 0 0 -72 -96 0 0 

0 0 108 0 -54 72 -54 -72 

k 0 0 0 192 72 -96 -72 -96 

150 _54_72-54 72 233 0 -125 0 

-72-96 72 -96 0 192 0 0 

0 0 -54-72-125 0 179 72 

, 0 0 -72-96 0 0 72 96 




and the corresponding compatibility conditions, using the method described 
in Section 2.1, become: 


0 — F X1 + F X2 + F X3 + F X4 (33) 

0 = F in +F va + F m + F w (34) 

0 = 4/ 5F X1 — 3/ 5F yi (35) 

0 = 4/5 F X2 + 3/5 F y2 + 6/5 F y „ (36) 


where (33) and (34) express the equilibrium of horizontal and vertical forces, 
respectively, and (35) and (36) establish the equilibrium of moments with 
respect to node 3 of the left and right substructures, respectively. 


Assuming that the forces acting on the structure satisfy the compatibility 
conditions (33)-(36) and take values: F = 1/8 (9 12 0 —8 8 —8 —17 4^j , the 
system (32) leads to the general solution: 



where p 1 ,p 1 ,p 3 ,p 4 G IR. From a mathematical point of view this is an affine 
space, i.e., the sum of a given vector (the first one) plus a linear space of di¬ 
mension 4 (an arbitrary linear combination of 4 linearly independent vectors). 
The graphical interpretation of these four vectors is shown in Figures 10 (a-d), 
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i.e., a horizontal translation, a vertical translation and two independent rota¬ 
tions with respect to node 3. Note that eliminating the bar 1-2 allows the bar 
1-3 and the right substructure to rotate with respect to node 3 independently. 
That is why there are two rotation parameters p 3 and p 4 . Note that in the 
solution (29) corresponding to the structure with this bar, the rotation vector 
is the sum of the last two vectors in (37). 



Fig. 10. Illustration of the four basic vectors of the linear space component of the 
solution corresponding to Example 8. They correspond to a horizontal translation, 
a vertical translation and two independent rotations with respect to node 3. 

If we add the boundary conditions U\ = 0,rq = 0 and u 4 = 0, then from (37) 
we obtain the system of linear equations: 


ui = 49/ (8k) + pi + 4/5p 3 = 0 
vi =p2~ 3/5 p 3 = 0 

U 4 = pi = 0 


which solution is pi = 0, p 2 = — 147/(32fc) and p 3 = —245/(327c), but p 4 
becomes free. 


Considering the boundary conditions matrix 



/I 0 0 0 0 0 0 0\ 



B = 

01000000 

i^ooooooioj 

yy rank 

; 


2 m —1=7 


that implies that there is 1 rigid solid movements allowed under the assumption 
of small displacements without strains. 
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Then, the solution in this case becomes: 


( U\ ^ 


( 0 ^ 


( ° ^ 

Vi 


0 


0 

U 2 


-3/2 


4/5 

V2 

1 

' k 

-13/2 

3 

+ Pa 

3/5 

0 

v 3 


-147/32 


0 

W4 


0 


0 



^ —147/32 ) 


\ 6/5 / 


that implies that the actual boundary conditions do not prevent the rotation 
of the right substructure with respect node 3 as shown in Figure 11. 



Fig. 11. Illustration of how the boundary conditions u\ = 0, v\ = 0 and = 0 do not 
prevent the rotation of the right substructure with respect node 3 in a mechanism 
structure. 



Fig. 12. The hyperstatic structure analyzed in Example 9. 


Example 9 (An hyperstatic structure) Consider now the structure given 
in Figure 12 where an additional bar joining nodes 1 and 4 has been added. 
The system of equations resulting from assembling the global element stiffness 
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matrices in this case is: 


Ku 


1 608 144 -500 0 0 0 -108-144^ 

144 567 0 0 0 -375-144-192 

-500 0 608 -144-108 144 0 0 

A; 0 0 -144 567 144 -192 0 -375 

300 o 0 -108 144 608 -144-500 0 

0 -375 144 -192-144 567 0 0 

-108-144 0 0 -500 0 608 144 

v —144—192 0 -375 0 0 144 567 y 


and the compatibility conditions become: 


H 


(F \ 

- 1 XI 

U 2 


F 

1 y i 

u 3 


F 

1 X2 

W 4 


F 

1 V2 

U 5 


F 

± X3 



F 

1 V3 

U 7 


f X4 

W 


K Fy4 J 


(38) 


0 — F Xl + F X2 + F x 3 + F Xa (39) 

0 = F vl + F va + F ya + F w (40) 

0 = 4/5 F X1 + 4/5 F X2 + 3/5 F m + 3/5 F m , (41) 


which correspond to the equilibrium of horizontal forces, vertical forces and 
moments with respect to node 3. 


Assuming that the forces acting on the structure satisfy the compatibility 

conditions (39)-(41) and take values F = 1/4 (3 8 0 —4 4 —4 —7 o) , the 
following general solution is obtained: 


( u x \ 


( 353/960 \ 


fl\ 


(0\ 


/4/5\ 

V\ 


521/540 


0 


1 


0 

U 2 


11/40 


1 


0 


4/5 

V2 

u 3 

1 

“ k 

-343/540 

133/192 

+ Pi 

0 

1 

+ P2 

1 

0 

+ P3 

3/5 

0 

V:i 


0 


0 


1 


0 



0 


1 


0 


0 



l o J 


\°) 




\w) 


Pi, Pi, P3 £ m. 


( 42 ) 

From a mathematical point of view it is an affine linear space, i.e., the sum of 
a given vector (the first one) plus a linear space of dimension 3 (an arbitrary 
linear combination of 3 linearly independent vectors). 


Example 10 (Hyperstatic structure) Consider the particular truss struc¬ 
ture in Figure 13, subject to 3 external forces. Note that it is the same example 
as the initial one but adding a new boundary condition (u 4 = 0), thus we add 
a new equation to the system (25) but we have a new unknown variable, the 
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Fig. 13. Hyperstatic truss structure of Example 10. 


vertical support reaction on node 4 (F yi ). Applying the superposition principle 
we transform our problem as it is shown in Figure 14. 

The set off all possible solutions will be composed by (29) and the particular 
solution associated with the vector 


F = - 

4 


-9 -4 0 0 0 0 9 4 


(43) 


which is 


Ad 

V\ 

u 2 

V2 

U3 

V'i 

U 4 

V 4 


\ U4 J 


( 467/120 ^ 
501/160 
359/120 
-359/160 
39/10 
0 
0 
0 


L _y± 

k 


^ —287/60 ^ 
27/40 
-179/60 
27/40 
-9/5 
0 
0 
0 


+Pi 


0 

1 

0 

1 

0 

1 

V 0 / 


+P2 


1 

0 

1 

0 

1 

0 

V 1 / 


+P3 


f 4/5 ^ 
-3/5 
4/5 
3/5 
0 
0 
0 


v u / v u / vj \ 6/5 / 

(44) 

where F yA1 p \, p 2 , p-^ G 1R , that from a mathematical point of view is an affine 
linear space, i.e., the sum of a given vector (the first one) plus a linear space 
of dimension 4 (an arbitrary linear combination of 4 linearly independent 
vectors). The graphical interpretation of the last three vectors is shown in 
Figure 5, whereas the first one is the particular solution associated with the 
system of forces (43). 


If we add the boundary conditions U\ = 0, V\ = 0,W4 = 0 and v 4 = 0, then 
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Fig. 14. Illustration of the superposition principle applied to Example 10. 
from (38) we have: 


( — 467/(120A:) ^ 

v 1 - 501/(160*;) 

«4 


V 


Va 




\ 


-287/(60 k) 10 4/5 ^ 
27/(40fc) 0 1 -3/5 
0 10 0 
0 0 1 6/5 J 


(F \ 
Pi 
P2 


(45) 


which solution is F Va = 317/269, p\ = 0, p 2 = 


\Ps / 

11267/(4304fc) and p 3 = 


56335/(25824fc). Thus the solution is unique and equal to: 


H 


0 1 

v-i 


0 

U 2 


657/538 

V2 

1 

-11867/4304 


- k 

957/538 

V3 


-11267/4304 



0 

w 


V 0 J 
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5 Solving Systems of Inequalities 


In many engineering applications we find systems of linear inequalities of the 
form: 

Ax < b. (46) 

Before trying to solve a system of inequalities of this type it is interesting to 
check whether or not the system is compatible, that is, if it has some solutions. 
Even, in many cases we can be interested in obtaining the conditions for 
the system to be compatible in terms of the independent terms b\, b- 2 , ■ ■ ■, b m . 
These conditions always have a physical or engineering meaning that raises 
some light about the problem under consideration. Once the system is proven 
to have solution we solve the system and obtain all its possible solutions. 

To solve these two problems the T algorithm is used (see Jubete [7,8], Padberg 
[9], Castillo, Jubete, Pruneda and Solares [5], Castillo, Esquivel y Pruneda [6] 
and Castillo, Conejo, Pedregal, Garcia and Alguacil [4]), that gives the dual 
cone of a given cone and is the key tool to discuss the compatibility problem 
and solve the system of inequalities. 

Since the concepts of cone and dual cone are used, we start with their defini¬ 
tions. 

Definition 1 (Polyhedral convex cone) Let A be a matrix, and {a 1; ..., a m } 

be its column vectors. The set 

A^ = {x G IR" | x = 7iqa! + • • • + 7T m a m with ttj > 0; j — 1,..., m} 

of all nonnegative linear combinations of the column vectors of A is known 
as the polyhedral convex cone generated by ai,..., a m (its generators), and is 
denoted A n . 

A cone A n can be written as the sum of a linear space V p plus a pure (acute) 
cone Wjr, i.e., A w = V p + W,. 

In this paper we use the Greek letter n to refer to non-negative real numbers. 

Definition 2 (Nonpositive dual or polar cone) Let A^ be a cone in IR" 
with generators ai,..., a*,. The nonpositive dual of A n , denoted A p , is defined 
as the set 

A p = {u G IR" | A T u < 0} = {u G IR" | afu <0; i = l,...,k} 

that is, the set of all vectors such that their dot products by all vectors in A n 
are nonpositive. 
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Definition 3 (Polytope) Let A be a matrix, and (ai,..., a. m } be its column 
vectors. The set 

{ m 

x G !R n | x = Aia] + • • • + A m a m with Xj > 0; j — 1,..., m; ^ A* = 1 

i=l 

of all linear convex combinations of the column vectors of A is known as the 
polytope generated by &i,, a, m (its generators), and is denoted A\. 

5.1 Deciding whether or not a system of linear inequalities is compatible 

In this section we show how to analyze the compatibility of a system of linear 
inequalities. 

First, we discuss the compatibility of a particular system of the form: 

a±iXi + 012^2 + ■ ■ ■ AainXn = bi, 

®2]Al +022^2 T ■ ■ ■ Aa 2 n X n = 62, 

. ■■ ■ ( 47 ) 

(In>\-1'\ 5~a m 2X2 A ■ 5~a mn x n b m 
xi,x 2 , ■■■ ,x n > 0 

that can be written as 



x 1 ,x 2 ,-" ,x n >0 


Expression (48) shows that the given system is compatible if and only if the 
vector b = ( 61 ,, b m ) T belongs to the cone generated by the set of column 
vectors (ai, a 2 , • • • , a n } of the coefficient matrix A, i.e., 

b G A w = b G (A£) p (49) 

Thus, the compatibility problem reduces to finding the dual cone V p + of 
the cone generated by the columns of the coefficient matrix and checking that 
b T V = 0 and b T W < 0. 
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To analyze the compatibility of an arbitrary system of linear inequalities, it can 
be converted to the case in (47), using slack variables to convert the inequalities 
in equalities, and one more artificial variable to convert the arbitrary variables 
into no negative variables, that is, each variable x t can be converted to x* — Xq, 
where xq, x* > 0. 


Example 11 (Compatibility of a linear system of equations in restricted variables) 

Consider the following linear system: 


x 2 -X'.i ~ 2 x 4 = b 4 
x 3 +x 4 = b 2 
-x 2 +x 3 + 2 x 4 = b 3 
x 2 +x 3 —x 4 = b 4 
—x\ +2x 2 +x 3 +x 4 = b 5 
Xi,x 2 ,x 3 ,x 4 > 0. 


(50) 


For it to be compatible, the vector b = (bi,b 2 ,b 3 ,b 4 } b 5 ) T must belong to the 
cone Cjr generated by the columns of the coefficient matrix, that is to say, it 
must belong to the dual of the dual of C n . Thus, the compatibility problem 
reduces to finding the dual cone = \ p + and checking that h 1 V = 0 
and h 1 W < 0. 


Since C v n is the cone: 



we obtain the desired compatibility conditions: 



b - 2 b 3 b 4 b 


( —2 1 —1 — 4 \ 

-3 1 -2 -7 

0 0 0 0 
1-112 
v 0 0 0 \) 


< 0 , 


( 51 ) 



b 2 b 3 b 4 b 


10 10 0 


b\ + b 3 


0 . 
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5.2 Solving a homogeneous system of linear inequalities 

Consider the homogeneous system of inequalities 

Ax < 0 (52) 

which can be written as 

a*x T <0; i — 1,..., m. (53) 

Expression (53) shows that (aq,... ,x n ) is the dual cone of the row vectors 
{a 1 , a 1 ,..., a” 1 } of A. 

Thus, obtaining the solution of the system (52) reduces to determining the 
dual cone A£ of the cone generated by the cone generated by the rows of 
matrix A. 

Thus, the general solution of an homogeneous system of linear inequalities is 
a cone, that is, its general solution is of the form 

p q 

X = £ PiVi + Y T’jWji pi e D': 7 Tj G ]R + . 

i=i j=\ 

Example 12 (Solving an homogeneous system of linear inequalities) 

Consider the system of equations 

—x 5 < 0 

X\ +X 2 +X 4 —x$ < 0 

—xi +x 2 -2 x 3 -x 4 +x 5 < 0 (54) 

—2x 4 +x 3 —x 4 —x 5 < 0 

2xi +x 2 —x 3 +x 5 < 0 . 

To solve this system, we obtain the dual cone of the cone generated by the 
rows coefficients and obtain the solution: 

2 2 0 -2 -2 \ 1 
-5 -1 -1 -1 3 ^2 

-1 3-1-13 7 t 3 

-5-1133 

7T 4 

0 0 0 0 4 / 

\ ^5 
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5.3 Solving a complete system of linear inequalities 


Now consider the complete system of linear inequalities: 

Ax < b (55) 


where A e IR mxn and b e IR m 

Adding the artihcial variable x n+ i, the constraint x n+ \ = 1 and the redundant 
constraint x n+ \ > 0 (it is a key trick that allows the constraint x n+ i = 1 to 
be easily forced at the end of the process), it can be written as 

anX\ +a± 2 X 2 + • • • -\-ai n x n — b\x n +\ < 0 

0'2lXI -\-a 22 X 2 + • • • +U2n^n —^2^n+l ^ 0 


5~a m 2X2 A ~\~a mn x n b m x n + 1 ^ 0 (^®) 

^ml^l S~Q>m2X2 "f" ' ' S~a mn X n b m X n -|_i ^ 0 

-X n +1 < 0 
Xn +1 1- 


System (56) can be written as 


(q-ii ,..., Q-in, b 1 ) (aq,..., x n , x n _|_i) ^ 0 

(021 j • • ■ » ®2 m ^ 2 ) (a?l 1 ■ ■ ■ 1 Xnt 2^n+l) 0 


( 1 j • • ^m) (^1) • • • i X n , X n +\) ^ 0 


-X n+ i < 0 


*^n+l 1 ■ 


(57) 


Expression (57) shows that (aq,,.. ,x n ,x n+ i) belongs to the dual cone of the 
cone generated by the set of vectors 

{ (^11 5 • • • i 11 in- ^1)5 (^21 5 • • • 7 ^2 m ^2)? • • • ? (^ml? • • • 5 dmm ^m)5 (h, 0, ■ 1). 1) } • 


Then, it is clear that the solution of (56) is the intersection of that cone 
with the hyperplane x n+ \ = 1. Thus, the solution of (55) is the projection on 
Xi x • • • x X n of the solution of (56). In other words, the general solution of 


29 




a complete system of linear inequalities is a polyhedron, that is the sum of a 
linear space, a cone and a polytope, that is, its general solution is of the form 

p q r r 

x = Pi w i + pi e IR, 7 tj e 1R + ; \ k = 1- 

i=l j=l fc=l fc=l 

Example 13 (Solving a complete system of linear inequalities) To solve 
the following system of inequalities: 


X\ 

+x 2 

+X 4 < 

1 

-Xi 

+X 2 -2x 3 

—X 4 < 

-1 

— 2 x\ 

+X 3 

—Xi < 

1 

2 xi 

+x 2 —X 3 

< 

-1, 


(58) 


we use the auxiliary variable x$ and the redundant constraint 1 = x$ > 0. 
Then, the system (58) can be written as: 






-X 5 

< 

0 

X\ 

+x 2 


+Xi 

-x 5 

< 

0 

—X\ 

+x 2 

-2x 3 

—Xi 

+X 5 

< 

0 

— 2 xi 


+x 3 

—Xi 

-X 5 

< 

0 

2 xi 

+x 2 

-x 3 


+2^5 

< 

0 





x 5 

= 

1 . 


(59) 


Since the upper part is an homogeneous system, one need to find the dual 
cone of the cone generated by the row coefficients. After imposing condition 
X 5 = 1 one gets the solution: 


f Xl ) 


^ —1/2 ^ 


1 2 2 0 -2^ 



x 2 

_ 

3/4 

+ 

-5 -1 -1 -1 


vr 2 

X 3 


3/4 


-1 3 -1 -1 


7T 3 



l 3/4 J 


v -5-l 1 3 ) 




which has no linear space part and whose polytope part reduce to a single 
vector. 
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6 The Truss Analysis Problem Revisited 


In this section we analyze the truss analysis problem, but assuming additional 
constraints related to node displacements and maximum bar compressions. 



Fig. 15. Truss structure of Example 14 showing the acting forces and the support 
conditions. 

Example 14 (The lifting jack structure) Consider the structure shown 
in Figure 15 where there is a lifting jack that push the node 4 in the vertical 
direction. We want to obtain the minimum force that has to be applied by the 
lifting jack, for the descent of node 4 to be smaller than 5 /k. We would like 
to get the set of all possible node displacements for all the valid values of the 
force exerted by the lifting jack. 

As the boundary conditions in this case are U\ = 0, V\ = 0, w 4 = 0 and 
v,i > —5//c, from (44) we can get the following system of linear inequalities: 

Ul - 467/(120 k) = -287/(60 k)F y4 + Pl + 4/5p 3 
vi - 501/(160 k) = 27/(40 k)F yi + p 2 - 3/5 p 3 

(61) 

W 4 = Pl 

5 /k > -p 2 - 6/5 p 3 

which solution, using the T-algorithm is: 


f f \ 

1 V4 


1 551/807 N 


( 1 ^ 

Pl 

P2 

= 

0 

-52441/(12912fc) 

+ 7r 

0 

233 / (807c) 

\ P 3 ) 


v -60595/(77472A;) y 


v 287/(48fc) y 


which is composed of a particular solution and a cone generated by a sin¬ 
gle vector. Note that F 1J4 = 551/807 is the force necessary to get a vertical 
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displacement of node 4 equal to —5 jk. Since it cannot be smaller and must 
satisfy the constraint u 4 > —5 /k, the value of n has to be positive. 

Replacing (62) in (44) we get: 


Cd 


0 ) 


0 ) 

Vl 


0 


0 

«2 


177/538 


9/5 

V2 

1 

-81521/12912 

7 r 

-l— 

287/40 

u 3 

“ k 

1437/538 

+ k 

-9/5 

V 3 


-52441/12912 


233/80 

U4 


0 


0 



v - 5 y 


v 807/80 y 


which, as before, is the sum of a particular solution and a cone generated by 
a single vector. 


From an engineering point of view, this solution must be interpreted as follows: 


(1) The particular solution gives the displacements produced by the force 
F yA = 551/807 that is the minimum force required to satisfy the con¬ 
straint i >4 > —5 jk. Note that the vertical displacement of node 4 is equal 
to —5 jk. 

(2) The cone generator represents the displacement increments produced by 
a unit force acting on the lifting jack. Note that the vertical displacement 
at node 4 is equal to 807/(80A;), and that because the n -value has to be 
positive the vertical displacement of node 4 is always greater than —5 jk. 

(3) Note that both solutions satisfy the boundary conditions u\ — 0, v\ — 0, 
m 4 = 0. 


Example 15 (Maximum compression constraint) Consider the same struc¬ 
ture as in Figure 7. The aim is to find all possible values of forces Pi, P 2 , and P 3 , 
in such a way that bars 1 and 5 are always subject to positive stresses (trac¬ 
tions). To solve this problem a relation between the displacements uf\vf\ 
rip , vj } and the axial force N ® acting on bar l must be established. From 
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(15) and (20) we have: 


N (l) = k (—1 0 1 o)l® u W = k 


cos cc 


(0 


„(0 


sina w cos a 


(0 


sm a 


(0 


(0 
'll' ' 


u 


(0 


) 

(64) 


If bars 1 and 5 are to be subject to tensions (Ad 1 ) > 0 and Ad 5 ) > 0) and the 
structure must satisfy the boundary conditions u\ = 0 , tq = 0 and w 4 = 0 , the 
following system of inequations obtained from (30) and (64) must be solved: 


Ui = 3Pi/(5fc) + 9P 2 /(20P) + 341P 3 /(120P) + 287P 4 /(60A;) + Pl + 4/5p 3 
Vl = 9Pi/(20k) + 19P 2 /(10P) + 25P 3 /(32/c) - 27P 4 /(40P) + p 2 - 3/5p 3 
U 4 = Pl 

o > 3P 3 /4 + 3P 4 /2 
0 > Pi + 3P 2 /4 + 3P 3 /2 + 3P 4 /2, 

(65) 

where the last two inequalities represent the constraints —Ad 1 ) < 0 and 
—Ad 5 ) < 0, respectively. 

The solution of the system using the T-algorithm is: 




^3/4^ 

p 2 


1 

p 3 


-2 

p 4 

= P4 

1 

Pl 


0 

P 2 


0 



1° J 


+P5 


/ 12 \ 


( 179 \ 

25 


f 1 l 


125 

16 




72 

_ 25 


4/3 


125 

0 


-4/3 


0 

0 

+7T 4 

0 

+7T 2 

0 



0 



0 




0 

1 


0 



~k 


233 


0 

rv 

, 0 ; 


V(72 k)J 


3 


; p 4 , Ps G 1R ; 


7Ti, 7T 2 > 0 


(66) 
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Replacing (66) in (30) we get: 




V\ 

U 2 

V'2 

U 3 

V 3 

W 4 

\ 1 ’ 4 ) 


P± 

k 


0 ^ 

0 

0 

25 

16 

0 

0 

0 

V 0 ) 


+t 


0 

0 

1 

0 

1 

0 

V 1 / 


+ T 


< o \ 
0 
1 

179 

~6CT 

0 

0 

0 

233 

^~ 60 ~^ 


+ T 


< o ^ 


0 

9 

T \ 


0 

0 
9 

V 10 ) 


:p 4 ,P5 e IR; 7Ti, 7T 2 > o, 


(67) 

which is a cone, i.e., the sum of a linear space of dimension 2 and an acute 
cone generated by two vectors. 


From an engineering point of view, this solution must be interpreted as follows: 

(1) The fist linear space generator in (67) corresponds to a vertical displace¬ 
ment of node 2 with the remaining nodes being fixed, which leads to 
compressions in bars 3 and 4 and no stresses in the remaining bars (see 
Figure 16(a)). 

(2) The second linear space generator in (67) corresponds to a vertical dis¬ 
placement of nodes 2, 3 and 4 with node 1 remaining fixed, which cor¬ 
responds to a rigid vertical displacement of the substructure defined by 
those nodes, while bar 2 is subject to tension stress. This implies that 
bars 1, 3, 4 and 5 are subject to no stress (see Figure 16(b)). 

(3) The first acute cone generator in (67) corresponds to a rotation with 
respect to node 3 of the substructure defined by nodes 2, 3 and 4, with 
the bar 2 remaining fixed, while bar 1 is subject to tension stress. This 
implies that bars 2, 3, 4 and 5 are subject to no stress (see Figure 16(c)). 

(4) The second acute cone generator in (67) corresponds to a vertical trans¬ 
lation of nodes 3 and 4 and a horizontal displacement of node 3. This 
implies that bar 2 is subject to compression, bar 5 to tension, and bars 
1, 3 and 4 are subject to no stress (see Figure 16(d)). 

(5) Note that all vertices satisfy the boundary conditions Ui = V\ = w 4 = 0. 

(6) Note also that the rotation, horizontal and vertical displacements asso¬ 
ciated with the acute cone generators must have the sign indicated in 
Figures 16(c) and 16(d) (7 t 4 , 712 £ IR + ) for the bars 1 and 5 being subject 
to positive stress (tension), while the vertical displacements associated 
with the linear space generators can have any direction (p 4 ,p 5 e IR) 
because bars 1 and 5 do not work under these displacements. 
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Fig. 16. Illustration of the different basic vectors that generate the linear space and 
the cone components of the solution of Example 15, where the tension bars 1 and 5 
have been distinguished from the remaining bars. 


Example 16 (Tension, compression and deflection constraints) Consider 
the same structure as in Figure 7, but now bars 3 and 4 have been replaced 
by cables, that can only stand tensions, and bars 1, 2 and 5 are built using 
concrete, material specially suitable to support compressions. 

The target of this example consists of obtaining the values of the forces 
P\, P‘ 2 , Pi and P 4 that can be applied on the structure in such a way that 
the cables work under tension and the concrete beams under compression, 
respectively. In addition, the maximum vertical deflection of node 4 is limited 
to V 2 > —1/k. 


If bars 3, and 4 are to be subject to tensions ( N ® > 0 and N ( 4> > 0), bars 
1,2 and 5 to compressions (Ad 1 ), Ad 2 ) < 0 and Ad 5 ) < 0) and the structure 
must satisfy the boundary conditions U\ = 0, Vi = 0 and u 4 = 0, the following 
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system of inequations obtained from (30) and (64) must be solved: 


u x = 3Pi/(60 k) + 9P 2 /(20/c) + 341P 3 /(120A;) + 287P 4 /(60A;) + Pl + 4/5p 3 
v x = 9Pi/(20fc) + 19P 2 /(lOfc) + 25P 3 /(32 /c) - 27P 4 /(40£;) + p 2 - 3/5p 3 
m 4 = Pi 

0 > —3P 3 /4 - 3P 4 /4 
0 > —5P 2 /4 - 5P 3 /4 - 5P 4 /4 
0 > —5P 3 /4 - 5P 4 /4 
0 > 5P 4 /4 

0 > -Pi - 3P 2 /4 - 3P 3 /2 - 3P 4 /2 
1/A: > 9Pi/20 + 27P 2 /80 + 233P 3 /160 + 27P 4 /40 - p 2 - 3/5p 3 , 

( 68 ) 

where the last six inequalities represent the constraints N ^ < 0 , N^ < 0 , 
—N^ < 0 , —N^ < 0 , i\d 5 ) < 0 and — v 2 < 1/A:, respectively. 


The solution of the system using the T-algorithm is: 



where 


= i; ^ > 0; * = 1,..., 5, 
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(70) 



and replacing (69) in (30) we get: 
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(71) 

, 5, which is the sum of a polytope 



Fig. 17. Illustration of the different basic vectors that generate the linear space and 
the cone components of the solution of Example 16. 

From an engineering point of view, this solution must be interpreted as follows: 

(1) The fist vertex in (71) corresponds to a vertical displacement of node 2 
with the remaining nodes being fixed, which leads to traction stresses in 
bars 3 and 4 and no stress in the remaining bars (see Figure 17(a)). 

(2) The second vertex in (71) corresponds to a vertical displacement of nodes 
2, 3 and 4 with node 1 remaining fixed, which corresponds to a rigid 
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vertical displacement of the substructure defined by those nodes, while 
node 2 is subject to compression stress. This implies that bars 1, 3, 4 and 
5 are to no stress (see Figure 17(b)). 

(3) The third vertex in (71) corresponds to a vertical displacement of bar 
4 and a rotation with respect to node 1 of the substructure defined by 
nodes 1, 2 and 3, which leads to a compression of bar 5. This implies that 
bars 1, 2, 3 and 4 are subject to no stress (see Figure 17(c)). 

(4) The fourth vertex in (71) corresponds to a rotation with respect to node 3 
of the substructure defined by nodes 2, 3 and 4, with the bar 2 remaining 
fixed, while bar 1 is subject to compression stress. This implies that bars 
2, 3, 4 and 5 are subject to no stress (see Figure 17(d)). 

(5) The fifth vertex in (71) corresponds to no displacement at all, which 
implies no stress in all bars (see Figure 17(e)). 

( 6 ) Note that all vertices satisfy the boundary conditions Ui = v\ = u 4 = 0. 

(7) Note also that the maximum allowed vertical displacement in node 2 
{v ‘2 > — 1/k) corresponds to all cases with A 5 = 0 , i.e., V 2 = —(Ai + A 2 + 
A 3 + A 4 ) = —1. 


7 Conclusions 


The following conclusions can be derived from this paper: 

(1) A full understanding of real problems stated as systems of linear equations 
or inequalities requires both the mathematical and the engineering points 
of view that complement each other. 

(2) The compatibility conditions must be interpreted from an engineering 
point of view, which help to identify errors, omissions or possible discrep¬ 
ancies between the mathematical model and the reality being modeled. 

(3) The mathematical structures of the general solutions, linear spaces, cones, 
polytopes and mixed combinations of these three structures have clear 
engineering interpretations that are closely related to the real problem 
being modeled. 

(4) The generators of the solution set, i.e., the linear space generators (basis), 
the cone generators, and the polytope generators (vertices) have clear 
interpretations from an engineering point of view, and contains a valuable 
information on the general solution of the problem and its properties. 
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